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e-MERLIN observations at 5 GHz of the GOODS-N region: 
pinpointing AGN cores in high redshift galaxies 



m ; 

o 

Oh- 



D. Guidetti *\ M. Bondi\ I. Prandoni\ R.J. Beswrf, T.W.B. Muxlow^, 
N. Wrigley^, I.R. Smail^ I. McHardy^ 

^ INAF - Istituto di Radioastronomia, via Gobetti 101, 1-40129 Bologna, Italy 

^MERLIN/VLBI National Facility, Jodrell Bank Observatory, 

University of Manchester, Macclesfield, Cheshire, SKI 1 9DL, UK, 

^Institute for Computational Cosmology, Durham University, Durham, DHl 3LE, UK, 

"^School of Physics and Astronomy, University of Southampton, Southampton S017 IBJ, UK. 



o 
u 

o 



> 

oo 

o 
m 

o 
m 



% 



Received 



ABSTRACT 

We present 5 GHz e-MERLIN observations of the GOODS-N region at sub-arcsec 
resolution (0.2-0.5 arcsec). These data form part of the early commissioning observa- 
tions for the e-MERLIN interferometer and a pilot for the e-MERLIN legacy program 
eMERGE. A total of 17 sources were detected with S/N > 3. These observations 
provide unique information on the radio source morphology at sub-arcsec scales. For 
twelve of these sources, deeper 1.4 GHz MERLIN-fVLA observations at the same 
spatial resolution are available, allowing radio spectral indices to be derived for ten 
sources on sub-arcsec angular scales. Via analysis of the spectral indices and radio 
morphologies, these sources have been identified as AGN cores in moderate-to-high 
redshift (1< z <4) galaxies. These results have provided AGN (or AGN candidate) 
classification for six previously unclassified sources and confirmed the AGN nature of 
the rest of the sample. Ultimately the eMERGE project will image the GOODS-N 
region at 1.4 and 5 GHz with higher resolution (about 50 mas at 5 GHz) and down to 
sub-/ijy sensitivities. The unique combination of sensitivity and spatial resolution will 
be exploited to study star formation and AGN activity in distant galaxies. 

Key words: galaxies: evolution - galaxies: active - galaxies: starburst - cosmology: 
observations - radio continuum: galaxies 



1 INTRODUCTION 

The apparent relationship between the mass of a super- 
massive black hole (SMBH) an d that of the pressu re- 
supported spheroid hosting it (e.g. Magor rian et al.lll998l ) is 
one of the most clear empirical indications of the potential 
importance of active galactic nuclei (AGN) and accretion- 
driven feedback mechanisms in galaxy formation and evolu- 
tion theories. Indeed, t he peak of QSO activity has been es - 
tablished at z ~ 2 fe.g. lSciimid^ll968l : lHasinger et al.ll2005h . 
at the same epoch when star formation activity in the Uni- 
verse was also at its peak. These processes often co-exist 
in galaxies, with observational evidence for the presence 
of an embedded AGN in ul traluminous starburst galaxies 
l| Alexander et al.ll2005l . 120081 ') and in 20-30% of ^ ~ 2 mas- 
sive star-forming galaxies in the Great Observatories Ori- 
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gins Deep Surveys (GOODS. ICiavalisco et all 120041 ') fields 
(|Daddi et alll2007l : lMullanev et al.ll2012l '). 

It is an essential if galaxy formation theories are to be 
tested that an assessment of the importance of embedded 
AGN activity in star forming galaxies be made at moderate 
to high redshifts. As such we need to identify the contri- 
bution of obscured accretion activity in bolometric surveys 
of the evolution of star formation. In other words, we need 
to identify and separately trace the total (obscured and un- 
obscured) star formation in the bulges of individual high- 
redshift galaxies, the related nuclear activity, and any star 
formation occurring on larger scales within a disk. Only by 
doing this will we be able to obtain a bolometrically com- 
plete census of star formation and of the growth of galaxies 
and their SMBHs. However, a complete census of both star 
formation and AGN galaxy activity, especially at high red- 
shifts, is challenging due to dust extinction and gas obscu- 
ration by circumnuclear material. 

Deep radio surveys provide a powerful, obscuration- 



D. Guidetti et al. 



independent tool for measuring both star formation and 
AGN activity in high-redshift galaxies, hence tracing the 
apparently simultaneous development of the stellar pop- 
ulations and the black-hole growth in the first massive 
galaxies. Indeed, multi-wavelength studies of deep radio 
fields show a com posite population of star-forming galax- 
ies and AGN (e.g. Prandoni et al.ll200ll: ICiliegi et al.( | 2005 : 



Afonso et al.ll2006: Smolcic et al.'2008: Se vmour et al.ll200S 



Mignano et al. 2008 : Padovani et al. 2009 1. with the former 
dominating at the lowest flux densities (S < 50 — 100 /iJy). 
Interestingly, about half of the AGNs probed in the deep- 
est radio fields are characterised by relatively low radio-to- 
optical flux density ratios (i? < C 100), as compared to classi- 
cal radio-loud AGNs (see e.g. I Padovani et al.ll2009l ). This is 
a direct evidence of radio-quiet AGNs being present at /ijy 
flux levels. Such a result supports the idea th at radio-quiet 
AGNs are not necessarily radio silent (see iKukula et al.l 
119981 ). and paradoxically opens the perspective of studying 
the entire AGN population, including the radio-quiet com- 
ponent, at radio wavelengths. 



1.1 Identifying AGN and Starburst Components 

Radio-quiet AGNs share many properties with star- 
forming galaxies: they have similar radio luminosities 
|--j^g22-24 -y^jj^-i-j ^^^ similar infrared/radio flux ratios. In 
addition radio-quiet AGNs are typically characterised by 
Seyfert-2-like optical spectra, which are often difficult to 
distinguish from tho se of star-forming galaxies (see e.g. 
I Prandoni et al.l |2009|) . This makes it difficult to separate 
radio-quiet AGNs from star-forming galaxies, even with 
multi-wavelength information, although the avai lability of 
IRAC colors and/or X-ray data can help fiPrandoni et al.l 
l2009l : lDel Moro et aLllJoTa ). 

One of the most direct ways to pinpoint embedded AGN 
cores in galaxies is the detection of a high surface brightness 
radio core through high resolution - VLBI-like - observa- 
tions. In order to reliably separate radio structures on several 
scales, from AGN core/inner jets (^Ikpc) to nuclear/disk 
starbursts (>lkpc) in high redshift {z > 1) galaxies, 
sub-/ijy imaging sensitivity togethe r with milli-arcsec 
resolution is clearly crucial (see e.g. iKewlev et al.l I2OO0I: 
Garredbool: iBiggs fc Ivisonll2008l : IChi. Barthel fc Garrett! 



20131 : iParra et al.ll2010l ). Pioneering Multi-Element Radio 



Linked Interferometer Network (MERLIN) observations 
at sub-arcsec resolution of the Hubble Deep field North 
(HDF-N) field, included in the GOODS-N region, did show 
that below a flux density of about 70 /iJy, the majority of 
the radio sources con sist of powerful star- forming galaxies, 
typically aX z < 1.5 (JMuxlow et al.l |2005| . hereafter M05). 
However a radio- AGN population was also morphologically 
identified, which had been classified as starbursts at all 
other wavelengths (jCasev et al.ll2009l ). 



1.2 e-MERLIN Pilot Observations 

The e-MERLIN interferometer is a major upgrade of the 
existing MERLIN array working at centimetre wavelengths. 
The array consists of up to seven radio radio telescopes 
spread across the UK and is equipped with a new optical 



fibre network and correlator, providing a ten-fold sensitivity 
boost, compared to MERLIN. 

The e-ME RLIN Galaxy Evolut ion Survey (eMERGE) 
Legacy project (jMuxlow et al.l 120081 ) represents the largest 
element of the recently approved e-MERLIN Legacy Pro- 
gramme, a set of large and high-impact surveys addressing 
fundamental questions in astronomy and astrophysics. The 
eMERGE Survey project aims to exploit the unique com- 
bination of e-MERLIN sensitivity (sub-/ijy rms noise) and 
spatial resolution (50 mas at 6 cm, comparable to that of the 
Hubble Space Telescope), to study star- formation and AGN 
activities in high redshift (1 < z < 4) galaxies, with special 
focus on possible co-existence and co-evolution of the two 
phenomena. One of two fields targeted by the project is the 
GOODS North (GOODS-N), the most intensively observed 
area on the northern sky, which includes the Hubble Deep 
Field North (HDF-N). 

Here we present early e-MERLIN commissioning data 
at 5 GHz of a 10 arcmin field within the GOODS-N region 
at sub-arsec resolution in the framework of the e-MERGE 
project. One of the aims of these pilot observations is to pin- 
point possible problems in the upgraded e-MERLIN system. 
The images are the first deep, sub-arcsecond e-MERLIN at 
5 GHz of such region. These 5 GHz e-MERLIN observations 
can be considered as complementary to those previously ob- 
tained at 1.4 GHz with MERLIN and presented by M05 who 
combined MERLIN and Very Large Array (VLA) images at 
1.4 GHz for 92 radio sources with Si. 4 > 4 iiij (extracte d 
from the original 1.4 GHz catalogue, see iRichardd |2000| ). 
and located within a region of 10 x 10 arcmin^ centred at 
R.A. 12'' 36™ 49^:4, Dec. +62° 12' 58'.'00 (J2000). These im- 
ages (M05), were not corrected for the primary beam atten- 
uation, have angular resolution in the range 0.2-0.5 arcsec 
and reach a la sensitivity of about 3/^Jy beam~^. The M05 
data have recently re-imaged over a large contiguous area 
and primary beam corrected (Wrigley et al., in prep) reach- 
ing a 1(T sensitivity of 3.5 /iJy beam~^ at the pointing centre. 

The new pilot C-band study presented here is ~5 times 
shallower than the previous L-band observations of M05. 
This, together with the lack of short spacings information, 
makes our pilot observations best suited to pinpoint rela- 
tively bright and compact sources associated to radio emit- 
ting AGNs, rather than probe faint and/or large-scale radio 
emission. Combining these pilot 5 GHz data with the new 
1.4 GHz images (Wrigley et al., in prep) we derived for the 
first time spectral indices at sub-arcsec resolution. 

The paper is organized as follows. Sects. [2]and[3]present 
the 5 GHz e-MERLIN commissioning observations, the data 
reduction and the source detection details. The spectral in- 
dices derived on sub-arcsec scale between the pre-existing 
1.4 GHz data and the new 5 GHz ones are discussed in 
Sect. O together with the diagnostics used to identify AGN 
activity in the observed galaxies. Individual source details 
are given in Sect. 15. Il where the classification of each source 
is discussed. Conclusions and future perspectives are sum- 
marized in Sect. [6l 



2 OBSERVATIONS AND DATA REDUCTION 

As part of the commissioning activities of e-MERLIN, 
in July and August 2011 we observed a 10 arcmin di- 
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ameter region within the GOODS-N field, centred at 
R.A. 12'' 36™ 4ffi0, Dec. +62° 12' 47'.'98 (J2000), for a total 
of 156 hours in the frequency range 4.5-5 GHz (hereafter re- 
ferred to as 5 GHz). This position differs by about 1 arcmin 
from the pointing centre used by M05, but was chosen to 
maximise S/N at known source locations. 

A sub-array of five of the seven e-MERLIN telescopes 
were used throughout these pilot observations, with the 76- 
metre diameter Lovell and the 32m Cambridge telescopes 
being unavailable at the time. The absence of the Cam- 
bridge antenna reduces the maximum spatial resolution of 
our 5 GHz observations, making it comparable to that of the 
previous 1.4 GHz MERLIN-I-VLA observations (see M05). 
Furthermore, the omission of the Lovell telescope gives rise 
to an approximately homogeneous array with 25m diameter 
elements, yielding a predictable primary beam. 

The 5 GHz data were taken in four adjacent 512 x 
0.25 MHz intermediate frequency channels (IFs) and ob- 
served both circular polarizations, with correlator outputs 
integrated every second. The observations were conducted 
with a cycle time of 1.5 min on the phase calibrator and 
8 minutes on the GOODS-N target yielding a total on-source 
integration time of 106 hours. 



2.1 Calibration 

The data were inspected, flagged, calibrated and imaged us- 
ing the NRAO AIPS package, following standard procedures 
for wide bandwidth observations. 

Visibilities that were clearly corrupted by instrumental 
phase errors, telescope errors, or radio frequency interfer- 
ences were spectrally removed (channel by channel) by using 
the tasks SPFLG and idled. 

All the calibrators and target data were then com- 
bined into a single multi-source file and averaged down to 
128 X 1 MHz channels, considerably reducing the size of 
the dataset. The absolute fiux-density scale was calibrated 
through observations of sources 3C 286 (primary flux cali- 
brator) and 1407-1-284 (OQ208, secondary point-source cal- 
ibrator) with flux density boot-strapped to 3C 286. Assum- 
ing a point-source model and a spectrum within the observ- 
ing band for OQ208 we established a fiu x density equal to 
2.45 Jy at the central observing frequency (jBaars et al.lll977l 
scale) . 

The flux calibrator 1407+284 was also used for bandpass 
calibration whereas phase correction of visibilities were cali- 
brated by observing the calibrator source 1241+603, located 
2° away from the pointing centre. The average flux density 
for 1241+603 was measured to be 0.216+0.007 Jy, a value 
consistent with the flux reported by the VLBA Calibrator 
List (NRAOfl (0.21 Jy at 5 GHz). 

The visibilities were then re-weighted to take into account 
the different sensitivities of the telescopes within the e- 
MERLIN array and finally time-averaged (4 s), in order to 
further reduce the data volume and speed up the imaging 
process. The effects of frequency and time averaging on the 
derived source fluxes are discussed in Sect. [S] 
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2.2 Imaging 

We produced total intensity images of all 92 radio sources 
presented in M05. These were obtained by averaging the 
four available IFs (task imagr) and restoring clean com- 
ponents with a circular beam of 0.5 arcsec diameter. Im- 
age dimensions measure 128x128 pixels and a cellsize 
equal to 0.0625 arcsec was chosen generating fields some 
8 arcsec across. Whenever a detection was found within 
one of the 0.5 arcsec resolution fields, a higher resolution 
(0.2 arcsec) map was re-imaged. The only exception is source 
J123725+621128, a relatively extended double radio source 
(see last panel of Fig. [T] and source details in Sect. 15. ip . 
where higher resolution imaging was not attempted in order 
to probe the extended emission visible at 1.4 GHz (M05). 
We also imaged 10 unresolved sour ces from t he or iginal 
1.4 GHz VLA catalog produced bv iRichardd (|2000l ') (see 
also lBiggs &: Ivisonll2006l : iMorrison et al.ll2010l l with arcsec- 
scale fiux density > 300 /iJy, located within a distance 
of 10 arcmin from our field centre at 0.2 arcsec resolution. 
These sources lie in a region outside the usable primary 
beam of MERLIN so were not included in M05. 

All the images were cleaned until the peak of the residu- 
als was about three times the rms noise [a). Cleaning deeper 
risks so-called "clean bias", tha t may affect deep imaging i n 
presence of poor UV coverage (|Whitelll997l : ICondonlll998l '). 
Frequency channels at the edges of the passband are gener- 
ally characterized by very low sensitivity, therefore the low- 
est rms noise level in the final images is obtained by exclud- 
ing these channels. This results in an effective bandwidth of 
120 MHz within each IF. 

The average off-source noise level, measured before cor- 
rection for the primary beam, in our naturally-weighted 
images at 0.2 arcsec and 0.5 arcsec resolution is ~ 15 and 
~ 20/iJybeam~^ respectively. From image analysis we ver- 
ified that the noise is close Gaussian, however, the values 
reported above are about 2 times higher than expected from 
thermal noise. The cause of sensitivity degradation in these 
early e-MERLIN commissioning 5 GHz observations have 
been identified with elements of the signal path introduc- 
ing noise in the data signal. Following these observations, 
additional first stage amplification has been introduced to 
mitigate this effect. 

The images were corrected for the primary beam at- 
tenuation, adopting the parameters used for the VLA an- 
tenna at 5 GHz. This is a reasonable assumption since all 
the antennas participating to our e-MERLIN observations 
are effectively 25-metre dishes. The primary beam correction 
results in rms noise which increases as a function of distance 
from the pointing phase centre, varying from approximately 
20-30 ^Jy within 200 arcsec of the pointing centre, to 40- 
50 /.iJy at 300 arcsec before rapidly growing to several hun- 
dred micro-Janskys. 



3 SOURCE DETECTION 

At a 3(T threshold, a total of 17 sour ces were detecte d (12 
from the M05 sample and 5 from the iRichardd (|2000|) cat- 
alogue). As expected, all these sources have relatively high 
measured fiux densities (> 80—100 ^Jy), and all but two are 
consistent with being unresolved. In contrast, all the M05 
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Table 1. Radio sources detected by e-MERLIN at 5 GHz. Col. 1: source name; Col. 2: distance from the pointing centre; Col. 3: local 
rms noise after the primary beam correction; Col. 4: Right Ascension and Declination of the peak; Col. 5: signal-to-noise ratio; Col. 6 & 
7: integrated flux densities at 5 GHz and 1.4 GHz, respectively; Col. 8: redshift. All parameters refer to 0.2arcsec resolution images, with 
the exception of source J123725-f 621128, for which the resolution is O.Sarcsec. In italic are the sources with uncertain flux density due 
to the large primary beam and smearing corrections. 
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sources were resolved at 1.4 GHz. This is easily explained by 
the lack of short baselines in our 5 GHz e-MERLIN obser- 
vations, together with a factor of ~ 5 worse sensitivity than 
for the 1.4 GHz observations (coupled with the fact that ex- 
tended synchrotron emission is typically characterized by a 
declining spectrum with frequency). 

The 5 sources from the I Richards! l|200CI ) catalogue are 
located at larger angular distances from the phase centre 
and are brighter than those in the M05 sample. 

The measured source parameters are presented in Ta- 
ble[T]which include the distance from the pointing centre, the 
local rms noise value, the peak position, the signal to noise 
ratio (S/N), the integrated flux density at 1.4 and 5 GHz. 
The source redshift is included, where available. For unre- 
solved sources, the integrated flux density is equal to the 
peak brightness, obtained through a quadratic polynomial 
fitting (aips verb maxfit), whereas for the two extended 
sources, the integrated spectral flux density is derived by 
integrating over a circular region with radius of 5 pixels cen- 
tred on the peak flux position. The peak brightness of the 
sources is corrected for the combined bandw i dth a nd time 
smearing effects following ISridle fc SchwabI l|l989 ) which 
predicts that sources within a distance of 300 arcsec from 
the pointing centre (imaged with a resolution of 0.2 arcsec) 
have intrinsic peak brightnesses reduced by up to 10%, while 
at ~ 400 arcsec the smearing reduces the peak flux den- 
sity by up to ~ 30%. These smearing corrections, in com- 
bination with primary beam uncertainties, yields highly un- 
certain flux density measurements for sources lying beyond 
400 arcsec from the pointing centre. In these cases we only 
provide a rough estimate of their flux densities which, be- 
ing all unresolved, are set equal to their peak flux density 
measurements. 



Contour images at 5 GHz of the 17 detected sources are 
shown in Fig. [l] The radio contours are overlaid on the pre- 
existing 1.4 GHz images obtained at the same resolution. 
As mentioned earlier , the 5 de t ected sources selected from 
the VLA catalog by I Richards! (|2000|) were not imaged by 
M05. Therefore, for such sources we do not have a sub-arcsec 
1.4 GHz image. 

In general the source peak positions at 1.4 and 5 GHz 
are coincident within the errors, suggesting a common origin 
of the two emissions. The two most notable exceptions are 
sources J123624-I-621743 and J123650-I-620844 where a sig- 
nificant offset is found between the 1.4 and 5 GHz emission 
centroids. These particular sources will be discussed further 
in Sect. (5] 

Figure [2] shows our source detection expectation based 
on the flux densities of the sources at 1.4 GHz used here as a 
proxy of those expected 5 GHz, under the assumption that 
the dominant component of the radio emission is a compact, 
flat spectrum (q = 0) AGN core. Such fluxes have to be 
compared to our 3(t detection threshold. It is clear that only 
the brightest sources in the M05 sample can be in principle 
detected. This is mainly a consequence of the fact that in our 
maps the achieved noise is, as already mentioned in Sect. [21 
a factor two higher than expected. 



4 PINPOINTING AGN CORES IN GOODS-N 
GALAXIES 

As discussed in Sect.[T] a clear separation between AGN and 
star-formation processes in galaxy samples can be very diffi- 
cult, especially at high redshifts, where gas/dust obscuration 
can be very important and composite systems (where both 
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Figure 1. Radio contours of the sources detected at 5 GHz overlaid to the 1.4 GHz emission obtained by the combined MERLIN- VLA 
data (Wrigley et al., in prep). The images with blank backgrounds do not have 1.4 GHz MERLIN- VLA data available. The restoring 
beam is 0.2arcsec for all sources. The contours are drawn at -2.5, 2.5, 3, 4, 5, 7, 10, 15, 20 times the local rms noise listed in Table[T] 



phenomena co-exist) are likely to be more frequent. One of 
the safest way to classify extragalactic radio sources as AGN 
or starburst is through optical spectroscopy. Unfortunatefy, 
these kind of observations can be extremely time consum- 
ing or not possible at all due to the faintness of the targets. 
Moreover even spectral classification is not always clear, es- 
pecially in presence of low luminosity Seyfert-2 AGNs, whose 
spectra are not always easily distinguished from those of 
star-burst galaxies. 



4.1 AGN/Starburst discrimination methodology 

A multi-wavelength approach is typically used when at- 
tempting to separate AGN from star-forming galaxies, pos- 
sibly taking advantage of all available diagnostic tools. Each 
diagnostic test, taken by itself, is usually not sufficient to se- 
curely classify galaxies, but when two or more criteria can be 
verified, the classification can be considered more reliable. 

The following diagnostics have been used by M05 to 
classify the 92 sources in their sample: 



(i) The radio-FIR correlation - this correlation is remark- 
ably tight for star-formin g galaxies and holds over a wide 



range of luminositie s (e.g. IdeJong et al 



■ 19851: ICondonlll992l : 
Garrettll20oi r How- 



iBeswick et al.ll2()0a ) and redshifts (e.g. 
ever the radio-FIR correlation is not useful to distinguish 
star-forming galaxies from radio-quiet AGNs (see discussion 
in Sect.[l]). 

(ii) The X-ray luminosity - low redshift star burst galax- 
ies are soft-band (< 2 keV) dominated sources (|Ptak et al.l 
1 19991 . [Alexander et al.ll2002l ) with X-ray luminosities Lx < 
10^^ W. X-ray luminosities Lx > 10^^ W, on the other hand. 



usually indicate the presence of an AGN IjRichards et al.l 
2007). However, exceptions exist, like e.g. low luminosity 
radio galaxies with X-ray lumino sities in the range Lx '^ 



10'' 



10^^ W (lEvans et al.ll2006l ). 



(iii) The radio morphology - jets/lobes emission is a clear 
indication of the presence of a radio AGN. A sub-arcsec-scale 
bright compact core can be associated to an AGN, although 
the possibility of a compact starburst can not be excluded 
at a resolution of 0.2arcsec (corresponding to ^^^ 1 — 1.5 kpc 
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Figure 1. (continued). Radio contours of the sources detected at 5 GHz overlaid to the 1.4 GHz emission obtained by the combined 
MERLIN- VLA data (Wrigley et al., in prep). The images with blank background do not have 1.4 GHz MERLIN- VLA data available. 
The restoring beam is 0.2arcsec for all sources, except for source J123725-f 621128, where it is O.Sarcsec. The contours are drawn at -2.5, 
2.5, 3, 4, 5, 7, 10, 15, 20 times the local rms noise listed in Table [T] 



at 1 < 2: < 



4 seelMuxlow et al.ll2005l: iRichards et al.ll2007l : 



iBateiat et aill2012l for Arp 220: iBondi et aLll2012l for Arp 
299). The optical appearance, if observations with enough 
resolution are available, can be usefully coupled to the radio 
morphology to discriminate between an AGN or starburst 
origin of the radio emission. 

(iv) The radio spectral index - radio emission with spec- 
tral index a < 0.4 {S ~ t^~°') is hkely to be powered by 
an AGN. On the other hand, steeper-spectrum emission can 
be produced by either radio jets/lobes (and therefore still 
AGN dominated) or by supernova remnants (starburst dom- 
inated). Inspecting the radio morphology with an adequate 
resolution can allow discrimination between these two pos- 
sibilities. 

By applying these diagnostic tools, a classification as 
AGN/AGN-candidate or SB/SB-candidate was proposed by 
M05 for most (~ 82%) of their 92 ra dio sources. Nine add i- 
tional sources were later classified bv lRichards et al.l (|2007|). 

As discussed in Sect, [l] one of the most direct ways to 



pinpoint embedded AGN cores in galaxies is the detection of 
a fiat-spectrum high surface brightness self-absorbed radio 
core through high resolution multi- frequency observations. 
In this paper, we take advantage of our pilot sub-arcsec 
5 GHz e-MERLIN observations, to better assess the pres- 
ence of AGN cores in GOODS-N galaxies. By combining our 
5 GHz maps with the recently re-CLEANed and primary beam 
corrected MERLIN-I-VLA data at 1.4 GHz (Wrigley et al. in 
prep), obtained at same resolution, we derived 1.4-5 GHz 
spectral indices for 10 of the 17 detected sources. Spectral 
index information and radio morphologies were both used 
to confirm previous AGN classification, when available, and 
to provide a new one for the sources still unclassified. The 
spectral index information is particularly relevant here, first 
because the radio morphology alone does not necessarily al- 
low us to discriminate a compact starburst from an AGN 
core (see discussion above, item iii). Secondly, brightness 
temperature measurements would not help us in the source 
classification, as the spatial resolution of our observations is 
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not high enough to be sensitive to excess of 10^ K. Indeed, 
the angular resolutions of our maps (0.2-0.5 arcsec FWHM) 
give brightness temperatures in a range (700-5000 K). 

Our classification also takes advantage of any pre- 
existing multi-w avelength information , in line with those 
used by M05 and lRichards et all (|200it . 

Spectral index analysis of M05 (see item iv above) was 
based on 1.4 and 8.5 GHz radi o fluxes obtain ed from arcsec 
resolution images (3.5 arcsec. iRichardg|2000l ). whereas our 
spectral indices are instead derived between 1.4 and 5 GHz 
sub-arcsec resolution flux densities. The latter are therefore 
better suited to pinpoint flat-spectrum compact AGN cores. 

For the unresolved sources at 5 GHz, the spectral in- 
dex is derived using the peak flux density of the corre- 
sponding 1.4 GHz component, when the peaks at the two 
frequencies overlap. For the two sources (J123624-f621743 
and J123650+620844) with a significant positional offset 
between the peaks at the two frequencies (see discus- 
sion in Sect. [3| we do not provide a spectral index. For 
the two extended sources at 5 GHz (J123644-I-621133 and 
J123646-I-621404) we derived the total flux at the lower fre- 
quency integrating over the same region used to derive the 
flux density at 5 GHz. It is worth noting that for source 
J123725-I-621128 we detect only the emission from a region 
eastwards the core, and therefore the spectral index refers 
to this area. 

One caveat must be stated clearly: whereas the 1.4 and 
8.5 GHz observations were made relatively close in time (1- 
2 years) , our e-MERLIN 5 GHz observations are separated 
from the 1.4 GHz data by more than 10 years. Variability in 
the radio domain may therefore be an issue and should be 
considered in the interpretation of our spectral indices. 



5 RESULTS 

The derived spectral indices are presented in Table [2] to- 
gether with the previous det ermination bet ween 1.4 and 
8.5 GHz at arcsec resolution iJRicharda l2000l ) . For sources 
not previously imaged at sub-arcsec scale by M05, we derive 
a spectral index using the arcsec-scale 1.4 G Hz (and 8.5 GHz , 
where available) flux density measurements ()Richardsll200d ) 
and our sub-arcsec 5 GHz flux densities. Such spectral in- 
dices have to be considered as tentative, since sub-arcsec 
5 GHz flux densities may be underestimated with respect 
to those obtained at lower resolution. For reasons of con- 
servatism, such spectral index values are not included in 
Tabled 

In general there is a reasonable agreement between the 
two spectral indices, although we notice a tendency for our 
1.4-5 GHz spectral indices to be flatter than those derived 
between 1.4 and 8.5 GHz. This is expected, as sub-arcsec 
resolution spectral indices are likely to better probe self- 
absorbed AGN cores than those derived at arcsec resolution. 
This supports multi-frequency sub-arcsec radio imaging as a 
promising technique to pinpoint AGN cores in high-redshift 
galaxies. 

Table [2] summarises the pre-existing source classifi- 
cations ba sed on sub-arc se c rad io imaging obtained by 
M05 and [Richards et al.l ((20071) an d, where available , 
on higher resolution VLBI imaging IjGarrett et al.l l200ll : 
IChi. Barthel fc Garretti 120131 ). The additional five sources 



Table 2. Spectral index and source classification comparison. 
Col. 1: source name. Col. 2: previous radio classification based on 
sub-arcsec radio imaging. Col. 3: spectral index between 1.4 and 
8.5 GHz ((Richards 2003). Col. 4: spectral index between 1.4 and 
5 GHz. Col. 5: new classification. The 5 sources listed at the bot- 
tom do not belong to the M05 sample and do not possess a clas- 
sification based on sub-arcsec scale imaging. 



Source 


class'" 


<i 


C'li 


class'' 


J 123606-1-620951 


AGN 


> 0.56 


0.53±0.23 


AGN 


J123617-(-621540 


unci 


> 0.55 


0.19±0.22 


AGN? 


J123623-I-621642 


AGN 


0.63±0.07 


1.02±0.23 


AGN 


J 123624-1-62 1743 


unci 






- 


J123642-I-621331 


AGN-I-SB? 


0.94±0.06 


0.86±0.18 


AGN-I-SB? 


J123644-I-621133 


AGN 


0.30±0.05 


0.24±0.11 


AGN 


J123646+621404 


AGN 


-0.04±0.06 


-0.15±0.20 


AGN 


J 123650-1-620844 


SB 


> 0.80 




- 


J123652-I-621444 


AGN 


0.12±0.07 


-0.19±0.25 


AGN 


J123714-I-620823 


AGN 


0.15±0.08 


0.04±0.09 


AGN 


J123721-I-621129 


AGN 


-0.28±0.06 


-0.35±0.22 


AGN 


J123725-f621128 


AGN 


1.35±0.06 


0.90±0.22 


AGN 


J123640-I-622037 


AGN 






AGN 


J 123649-1-620438 


- 






AGN 


J123649-f620737 


- 


0.56±0.07 




AGN? 


J123659-I-621832 


- 


0.26±0.07 




AGN 


J123721-(-620708 


- 






AGN? 


" Garrett et al. (2001). M05, 


Donley et al. 


(2005), 




Richards et al. f2007~). Chi, Barthel & Garrett f2013): ' 


' This work. 



not belonging to the M05 sample (listed at the bottom) 
do not therefore have such classifications available. As ex- 
pected, eight of our sources are associated wi th galaxies clas- 
sified as AGN or AGN candidates by M05 or [Richards et al.l 
(2007) (Fig. [2ll. One of the additional sources was classi- 
fied as a possible composite starburst-|-AGN object (Fig. [2]). 
Six sources were also detected by higher resolution Very 
Long Baseline Inter f erome t ry (VLBI) observations at L - 
band (JGarrett et all l200ll : IChi. Barthel fc Garretd l2013l l. 
confirming the presence of an AGN core at the very cen- 
tre of the galaxy. Only one source was previously classified 
as a starburst galaxy. Our new classification is discussed for 
each source in Sect. ISTTl and the results are summarized in 
the last column of Table [2] 

We detect sources preferentially classified as 
AGN/ AGN- candidates rather than starburst compo- 
nents which, as mentioned earlier (see Sect. [3]), is a 
consequence of the fact that our 5 GHz observations are 
sensitive only to the bright compact radio sources in the 
observed field, and of the fact that AGNs are detected in 
large fractions at the brigh ter end of sub-mjy rad io samples 
(78% at S'i.4GHz > 500 uJv. lMignano et all (|2008D . Indeed of 
the 19 sources in our sample S'i.4ghz > 300 /xjy, 10 of the 11 
detected are positively identified as AGN/AGN candidates 
(Fig. [2]). Under the conservative assumption that the unde- 
tected sources are all lower surface brightness star-forming 
galaxies, we can conclude that at such flux densities ~ 53% 
(10/19) of the sources are AGN-dominated. Going to lower 
flux densities (100 < 5i.4ghz < 300 /xJy), and limiting our 
statistical analysis to sources within d < 300 arcsec from 
the image centre, where our detection threshold is sensitive 
enough (see Fig. [2]), we detect only 4 additional sources 
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Figure 2. 1.4 GHz flux density versus distance from the pointin g 
centre for the 92 sources imaged by M05 and iRicharda l|200CI ) . 
Undetected sources are shown as black empty circles for the M05 
sample and red filled circles for lRichardd l|2000|) sources. Detected 
sources are indicated by symbols in different col ors correspondin g 
to different source classifications (from M05 or iRichardj ||200CI ). 
see discussion in Sect. O: AGN (green triangles); starburst (SB, 
blue pentagons); possible composite source (AGN+SB?, orange 
pointed cross es); uncla s sified source (black circles and red stars 
for M05 and IRicharda l|200Cl ) sources, respectively). The black 
horizontal line indicates the flux density limit of the M05 sample 
(92 sources with S(i.4GHz) > 40/xJy), whereas the red horizontal 
line indi cates the flux de nsity limit of the additional 10 sources 
from the IRicharda l|200(J ) catalogue imaged in the framework of 
this work (Sn 4qhz) > SOOfiJy). The dashed and dotted lines 
represent the rms measured in the 5 GHz images we produced 
at 0.2 and 0.5 arcsec resolution, respectively. The dot-dashed line 
represents our Scr detection threshold. 



(out of 19 in these flux and distance intervals). Three 
of such sources are associated to AGN/AGN-candidates, 
corresponding to ~ 16% of the total (3/19). This has 
to be considered as a conservative estimate of the AGN 
fraction, as we cannot exclude, in this flux interval, to 
have steep-spectrum AGN among the undetected sources. 
Nevertheless our results are consistent with a decreasing 
fraction of AGNs with decreasing flux densi ty, as found in 
wider arcsec-scale sub-mjy samples (see e.g. ISevmour et al.l 
l2008l ). 



5.1 Source classification 

In this subsection we present the updated classifications 
based on our 5 GHz sub-arcsec imaging and spectral index 
analysis of each source. Each classification is discussed with 
respect to previous indications from M05 and other multi- 
wavelength observations. 

J123606+620951: This sour ce was unclassi f ied b y 
M05 and later classified as AGN by [Richards et all |20o3), 
thanks to new X-ray band information. The integrated flux 
density measured directly from the 0.2" resolution 1.4 GHz 
images is higher than that reported in Table A2 by M05 and 



derived from VLA only observations. It is possible that the 
source is variable. The core position at 5 GHz agrees within 
the errors with that were reported by M05 at 1.4 GHz. The 
spectral index we obtain is consistent with the arcsec-scale 
previous determination between 1.4 and 8.5 GHz, and is not 
particularly informative (a ~ 0.5 ± 0.2). However the de- 
tection of a compact, possibly variable, component likely 
confirms the AGN origin of the radio emission. 

J123617+621540: M05 observed a compact struc- 
ture together with two-sided emission extended over about 

1.3 arcsec, but the source was unclassified because the ra- 
dio spectrum was poorly constrained (a > 0.55). How- 
ever, we clearly detect a bright core with flat spectral index 
(a ~ 0.2±0.2), strongly suggesting the presence of an AGN. 

J123623+621642: The source was classified in M05 as 
an AGN candidate. The presence of an AGN core was later 
confirmed by global VLBI o bservations at 1.4 GHz, that re - 
vealed a compact emission (|Chi. Barthel fc Garretj 120131 1. 
At 5 GHz we detect a compact component coincident with 
the peak of the 1.4 GHz radio source. The spectral index we 
derive is steeper (a ~ 1.0 ± 0.2) than that found between 

1.4 and 8.5 GHz, suggesting possible flux density variabil- 
ity between the 1996-97 1.4 and 8.5 GHz observations and 
our 2011 5 GHz observations. The compact structure and 
the possible variability both confirm the AGN origin of the 
radio emission. 

J123624+621743: The source does not have a 1.4 - 

8.5 GHz spectral index measurement. It has a faint X-ray 
counterpart but the lack of redshift information does not 
allow a derivation of the X-ray luminosity. We detect a weak, 
compact structure at 0.2 arcsec resolution with S/N of 3.3, 
that is not coincident with the peak of the 1.4 GHz image, 
but overlaps to its northern extended radio emission. To 
check the reliability of this detection, we used the Gaussian 
statistics to evaluate the probability to have a noise spike 
in the 5 GHz image above 3.3a within the area covered by 
the 1.4 GHz VLA emission. This assumption is motivated 
by the fact that the rms noise distribution is very close to 
Gaussian. The probability for a spurious signal is of about 
0.7%, however deeper observations are needed to confirm 
this detection. 

Because of this offset, we have not calculated the spec- 
tral index. The source is associated with a relatively faint 
1.4 GHz source (5(1.4 ghz) ~ 70/iJy) and was not expected 
to be detected at 5 GHz (it is well below our 5 GHz detection 
threshold in Fig. [2]). The presence of a significant positional 
offset between the two peaks at 1.4 and 5 GHz and the appar- 
ent flux inconsistency make this detection uncertain. Deeper 
observations arc therefore needed to confirm this detection. 

J123640+62203T: This is one of the sources not pre- 
viously imaged at 1.4 GHz by M05, and therefore a classi- 
fication was not previously provided. However this object 
shows a radio-excess with respect to the radio/infrared cor- 
relation for radio-quiet AGNs and star-formin g galaxies and 
was therefore classified as a radio-loud AGN l|Donlev et al.l 
|2005|). This source is neither detected by Chandra nor has 
a 1.4 — 8.5 GHz spectral index determination. At 5 GHz we 
observe a bright compact core with a highly uncertain flux 
density due to the large primary beam correction. With the 
caveat of this uncertainty, we tentatively obtain a radio spec- 
tral index between 1.4 and 5 GHz of ~ 0.4. The spectral in- 
dex may be flatter, if we consider that some of the 5 GHz 
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flux could be resolved out in our sub-arcsec 5 GHz images. 
This seems to confirm the AGN classification inferred by the 
high radio/infrared flux r atio. 

J123642+621331: IWaddington et al.l (|l999l ') argued 
that this source is most likely a distant {z = 4.424) 
dusty starburst with a weak embedded AGN. At 1.4 GHz 
the source shows both diffuse and compact structures 
(M05). The latter has been reveale d at higher resolu- 
tion by VLBI observations at 1.4 GHz (iGarrett et al.ll200ll : 
IChi. Barthel fc Garretti |2013| ) , confirming the presence of 
AGN emission. At 5 GHz, we detect a compact component 
coincident with the peak of the 1.4 GHz emission. Since at 
1.4 GHz the radio source is extended, we use only the flux of 
the compact emission at 1.4 GHz to derive a spectral index 
of a ~ 0.9 ± 0.2. This value is consistent with the spectral 
index derived at arcsec-scale between 1.4 and 8.5 GHz. The 
steepness of the spectrum may conflrm the composite nature 
of the source, and/or be explained by steepening effects in- 
troduced by the (large) fc-correction of this high redshift 
source. 

J123644+621133: This is a FR I radio galaxy. 1.4 GHz 
sub-arcsec scale imaging shows a bright compact core and a 
short jet extending southwards (M05). 1.4 GHz VLBI obser - 
vations (JGarrett et all 1200 ll : IChi. Barthel fc Garrett 120131 ') 
detected the compact core component. A core-l-jet morphol- 
ogy is observed also in our 5 GHz image, conflrming the AGN 
nature of the source. Considering the same area emitting 
at 5 and 1.4 GHz, we derive a rather flat spectrum with 
Q~ 0.2 ±0.1. 

J123646+621404: An unresolve d core has been de- 
tected by VLBI observations at 1.4 GHz (JGarrett et al.ll200ll : 
IChi. Barthel fc Garrettll2013l l. Extended emission southeast 
of the nucleus was detected at 1.4 GHz by M05. At 5 GHz 
we observe a compact core and again some extended emis- 
sion southeast of the nucleus. The radio spectrum is inverted 
on both arcsec and sub-arcsec scales, conflrming the AGN 
nature of the radio emission. 

J123649+620438: This is one of the sources not pre- 
viously imaged at 1.4 GHz by M05 which does not have a 
1.4 - 8.5 GHz spectral index determination. The AGN/SB 
discriminators in this source are ambiguous. The galaxy is 
detected by Chandra with a 0.5-8 keV X-ray luminosity of 
Lx — 5.3 X 10^'^ W, and a rather soft-band dominated emis- 
sion, r = 1.6, that wo uld be consistent with starburst emis- 
sion ([Alexander et al.ll2002 V On the other hand, on the basis 
of an optical spectrum showing the Ha line, the g alaxy was 
classifled as an AGN (|Hornschemeier et al.l |2005|) . We de- 
tect a bright core at 5 GHz, that implies a flat or slightly 
inverted radio spectrum between 1.4 and 5 GHz (q ~ —0.1), 
supporting the low luminosity AGN scenario. 

J123649+620737: This is one of the sources not previ- 
ously imaged at 1.4 GHz by M05. The overall radio spectral 
index on arcsec-scale is steep {a ~ 0.56 ± 0.07, see Table [2]), 
and the 1.4 GHz arcsec-scale image shows extended emis- 
sion. The radio source is associated with a hot-dust domi- 
nated ultraluminous galaxy thou ght to host both a n AGN 
and intense star forming regions (jCasev et al.ll2009l ). Alter- 
natively, the extended 1.4 GHz emission may be AGN re- 
lated extended emission and the 70 nm emission may be 
due to a dusty torus surrounding the AGN. The compact 
component detected at 5 GHz overlaps with the compact 
bright optical nucleus visible in the HST image (see Figure 



3 in lCasev et al.ll2009l V Using our 5 GHz flux density mea- 
surement we flnd a possible steepening of the radio spectral 
index: a ~ 0.5 between 1.4 and 5 GHz and a ~ 0.8 between 
5 and 8.5 GHz. If conflrmed, this is consistent with optically- 
thin synchrotron emission, associated with the central AGN. 

J123650+620844: This source is classifled as a star- 
burst by M05. It does not have a 1.4 — 8.5 GHz spectral index 
determination. A weak unresolved core is detected at 5 GHz 
with S/N of 3.5 The peak is slightly offset (0.18±0.04 arcsec 
northwards) with respect to the 1.4 GHz peak position, but 
still overlaps to some extended emission detected at the 
lower frequency. As with J123624-)-621128 we have not cal- 
culated the spectral index. The presence of an offset may 
support the starburst classiflcation. As for J123624-f 621743, 
we estimated the probability to have a noise spike within the 
area corresponding to the 1.4 GHz VLA emission above 3.5a. 
In this case, the probability for a spurious signal is of about 
0.3%. However, future observations with higher sensitivity 
and better M^^coverage are needed to confirm this detection. 

However, as for J123624-I-621743, the source was not 
expected to be detected at 5 GHz: this source has a flux 
S'(i.4GHz) ~76/iJy and is located well below our 5 GHz de- 
tection threshold (see Fig. [2]). The presence of a signiflcant 
positional offset between the peaks at 1.4 and 5 GHz and the 
apparent flux inconsistency make this detection doubtful. 

J123652+621444: This source is characterized by a 
flat arcsec-scale spectral index and was classifled as an AGN 
by M05. The radio structure at 5 GHz shows a compact core 
coincident with the peak at 1.4 GHz. The derived spectral 
index is inverted (q ~ —0.2 ± 0.2), conflrming the AGN 
nature of the source. 

J123659+621832: This is one of the sources not pre- 
viously imaged at 1.4 GHz by M05. The radio source is 
identifled with a dust-obscured galaxy, an objects show- 
ing an IR excess with resp ect to the optical emission 
(IGeorgantopoulos et al.l |201l|). The object is detected by 
Chandra and shows evidence for moderate absorption. 
The redshift is not constrained: a tentative estimate of a 
phot ometric redshift yields a value of 4.4 (jRovilos et al.l 
|201Q^: in any case the sour ce must be at redshift > 2 
(IGeorgantopoulos et al.ll201lh . The arcsec-scale spectral in- 
dex is flat. At 5 GHz we detect a very bright compact source 
which implies an inverted radio spectrum between 1.4 and 
5 GHz (a ~ —0.4), and a very steep spectrum between 5 and 
8.5 GHz (a ^ 1.9). Such spectral behaviour might be due to 
variability. However this kind of spectrum i s also typica l of 
Giga-Hertz Peaked radio sources (GPS, e.g. lQ'Dealll998h . a 
class of very compact radio galaxies, thought to be caught 
in the very early phases of their growth. Whatever the case, 
this points towards an AGN classification. Assuming a con- 
servative redshift of z ~ 2 and a conservative radio spec- 
tral index of a ~ 0.26 (the one measured between 1.4 and 
8.5 GHz at arcsec scale, see Table [2|, we derive a rest-frame 
frequency fc-corrected luminosity of 3 x lO'^^ WHz~^. Even 
under such conservative assumptions, this turns out to be 
the highest radio luminosity source in the sample. This po- 
tentially very high-redshift source is of particular interest 
and worth further investigation. 

J123714+620823: The core was detected at 14 GH z 
by global VLBI observations (|Chi. Barthel fc Garret3l2013l '). 
confirming the AGN classification by M05. The arcsec-scale 
radio spectrum is fiat, indicating that the overall radio emis- 
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sion is dominated by the AGN core. This is consistent with 
our own detection of a bright compact core, with a flat radio 
spectrum (a ~ 0.04 ± 0.09) between 1.4 and 5 GHz. 

J123721+620708: This is one of the sources not previ- 
ously imaged at 1.4 GHz by M05 and does not have a 1.4—8.5 
spectral index determination. Our 5 GHz detection on sub- 
arcsec scale implies an inverted radio spectrum of a ~ —0.5 
between 1.4 and 5 GHz. This may be due to variability, sup- 
porting an AGN classification. However the significance of 
the 5 GHz detection is low (3cr) and deeper observations are 
necessary to confirm the AGN nature of this source. 

J123721+621129: The core was detected at 14 GH z 
by global VLBI observations (|Chi. Barthel fc Garret3l2013l ). 
confirming the AGN classification by M05. The arcsec-scale 
radio spectrum is inverted, indicating that the overall radio 
emission is dominated by the AGN core. At 5 GHz we also 
detect an inverted-spectrum (a ~ —0.4 ± 0.2) unresolved 
core, that overlaps the peak emission at 1.4 GHz. 

J123725+621128: This is a wide-angle tailed (WAT) 
radio galaxy. 1.4 GHz sub-arcsec scale imaging shows a rel- 
atively extended double source. At 5 GHz we detected the 
eastern hot-spot/jet component of this source. The derived 
spectral index between 1.4 and 5 GHz is steep with a ~ 0.90 
and even steeper between 1.4 and 8.5 GHz (a = 1.36). The 
latter, being derived with a much lower resolution, likely 
contains significant contribution from the diffuse tail emis- 
sion. 



6 SUMMARY & FUTURE PERSPECTIVE 

The e-MERLlN commissioning observations presented here 
represent the first high resolution wide-field imaging of the 
GOODS-N field at 5 GHz. which. Although Umited to 5 an- 
tennas of the e-MERLIN array the observations yield a reso- 
lution of 0.2 arcsec and a depth of about 15 /iJy rms noise at 
the phase centre. These observations were used in conjunc- 
tion with archival data to determine classifications of 15 of 
17 sources detected within the field. 

We detected 12 out of the 92 sources belonging to 
the complete sample derived from deeper 1.4 GHz MER- 
LIN-I-VLA observations by M05. Comparing our 5 GHz im- 
ages with the MERLIN+VLA images at 1.4 GHz at the same 
resolution we derived for the first time sub-arcsec scale ra- 
dio spectral indices for 10 of the 12 sources detected. The 
remaining two sources have displaced peaks at the two fre- 
quencies and the 5 GHz detection needs confirmation. In ad- 
dition, we detected 5 more sources located outside the im- 
aged field in M05. For such sources our 5 GHz e-MERLIN 
observations provide unique information on the radio source 
morphology at sub-arcsec scales. 

Due to the still limited sensitivity of the present ob- 
servations and the lack of short baselines, our observations 
are sensitive to relatively bright compact radio emission, 
most likely of AGN origin, and indeed many of the detected 
sources have independent indicators of AGN activity. The 
analysis of 5 GHz morphologies, together with the derived 
spectral indices allowed us to confirm previous AGN clas- 
sifications, when available, and propose a classification as 
AGN/ AGN-candidate for six of the seven sources previously 
unclassified. Of eight radio sources previously classified as 
AGN, six have 1.4-5 GHz spectra flatter than 0.6. 



The fact that most of the detected sources are associ- 
ated with AGN/AGN-candidates confirms previous results 
showing that AGNs dominate at the brighter end of sub-mjy 
radio samples (e.g. Mignano et a l. (2008)). The relative frac- 
tion of detections (and AGN) decreases going to lower fiux 
densities (AGN sources go from 56% at Si.4ghz > 300 /^Jy, 
to > 16% at 100 < 5i.4GHz < 300 /xjy), where the compo- 
nent associated with steep-spectrum star-forming galaxies is 
expected to become more prominent. 

The detected sources span a very wide redshift range 
(0.1 < z < 4.4). Three of the highest-redshift sources 
in the detected sample have independent observational ev- 
idence (based on infrared/sub-mm data) that the AGN 
might be embedded in dusty star-form ing galaxies: they 
are J123642-I-621331 at 2 = 4.424 (|Waddington et al.l 
1999D. J123 649+620737 at z = 2.315 (jsmail et all |2004| : 
Casev et al.] [2009) and J123659-I-621832 at z > 2 
(ten tative photom etric redshift estimate of 2; ~ 4.4, 
.Rovilos et al.ll2010l ). The latter two are among the sources 
for which sub-arcsec radio imaging was not previously avail- 
able. Interestingly, the radio spectral behaviour of source 
J123659-I-621832 seems to mimic that of GPS radio sources. 
This potentially very high-redshift source is of particular in- 
terest and is worth further investigation. This pilot study 
provides a indication of what will be achieved in the near 
future now that e-MERLIN is fully operational. In the fol- 
lowing months, the GOODS-N region will be targeted at full 
resolution (about 50 mas at 5 GHz) by deep (sub-/xJy sensi- 
tivity) e-MERLIN 1.4 and 5 GHz observation s in the frame- 
work of the eMERGE Survey Legacy project (jMuxlow et al.l 
[200«). The 5 GHz observations will provide the resolution to 
directly distinguish between jet-like or more compact mor- 
phologies, indicating AGN activity, and less compact radio 
emission associated with star-forming galaxies. This analysis 
will benefit from 1.4 GHz information, allowing derivation of 
the source radio spectral index on sub-arcsec scales. More- 
over, ongoing complementary deep Karl G. Jansky Very 
Large Array (JVLA) observations will recover the extended 
emission unsampled by the e-MERLIN data on angular 
scales larger than 0.5 arcsec. The combination of deep, high 
resolution 1.4 and 5 GHz JVLA and e-MERLIN observations 
will provide the opportunity to both spectrally and morpho- 
logically identify embedded low-luminosity AGNs, whilst at 
the same time map the distribution of star formation activ- 
ity on scales from 50 — 2000 mas in the high-redshift galaxy 
population (1 < z < 5). These angular scales correspond to 
spatial resolutions spanning sub-kpc (~ 300 pc) to tens of 
kpc at z>l. This strategy will enable us to obtain a bolo- 
metrically complete census of star formation, the growth of 
galaxies and their SMBHs, an essential step toward a full 
understanding of galaxy formation and evolution. 
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